Extracting Phases from Aperiodic Signals 
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We demonstrate by means of a simple example that the arbitrariness of defining a phase from an 
aperiodic signal is not just an academic problem, but is more serious and fundamental. Decompo- 
sition of the signal into components with positive phase velocities is proposed as an old solution to 
this new problem. 
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There are some quantities which, although in principle 
exactly and uniquely defined, are hard to estimate from 
experimental data. Prominent examples are informa- 
tion theoretic quantities like algorithmic cornplexities (for 
which onlyupper bounds can be estimated 0]), Shannon 
entropies P| , and mutual entropies !3| . Opposed to these 
are quantities where even a unique definition is lacking, 
although most researchers believe that they have a decent 
common sense "definition" . Maybe the most important 
quantity in this category is the phase of a non-periodic 
signal nil HQ 

For a pure sine wave signal the notion of phase is 
obvious and trivial, provided the sampling rate is high 
enough (which we will assume throughout the following) . 
Things are nearly as clean for anharmonic periodic sig- 
nals. There, we can map the orbit, e.g. by delay embed- 
ding, onto a closed loop in a plane, and we can define 
the phase by the angle of the vector from some point 
in the interior of the loop to the point corresponding to 
the actual state. This phase will of course depend on 
the central point and on the delay, but different choices 
will give equivalent phases, if the loop does not intersect 
itself: They will give the same average angular veloc- 
ity uo = hmT^oo(0(i + T) — (f){t))/T, and the difference 
between two phases defined that way will stay bounded 
with time. This "geometric" definition of phase can be 
generalized to aperiodic signals whenever there is a way 

- either via embeddings or using multivariate time series 

- to project the orbit into a plane in such a way that 
it always encircles some point. A typical example is the 
Rossler attractor Q for particular values of its parame- 
ters. But if the loop intersects itself such that its interior 
is divided into several domains, then central points cho- 
sen in different regions will lead to non-equivalent phases 
(see Fig.l) 

In addition to this geometrical definition, another pop- 
ular approach is via the Hilbert transform 0,HI3- Under 
mild restrictions on the signal x{t), the pair {x{t),y{t) = 
{Hx){t)} form real and imaginary part of an analytic 
function, and the phase is defined as its argument in 
the complex plane. In general this gives a well-defined 
phase. Its value changes of course if the signal is shifted, 
x{t) — + x{t) + c, but a unique phase is obtained after 
de-meaning, i.e. when c is such that {x{t)) = after the 
shift. Problems arise if the trajectory goes through the 



FIG. 1: Projection of a periodic orbit into a plane such that 
it intersects itself and gives rise to 2 domains. For central 
points chosen in domain II the mean angular velocity is half 
of that for central points in domain I. 



origin, x{t) = {Hx)[t) = at some time t p|. Prac- 
tically the method breaks down even when this is not 
exactly fulfilled, since then the phase is very sensitive to 
low amplitude noise. 

A third situation where one might be tempted to see a 
"natural" way of defining a phase is when a signal arises 
from circular motion in some d-dimensional space with 
constant amplitude A and arbitrarily changing 4i{t), 



x{t) = Acos(0(t)) 



(1) 



Certainly this is considered by many as the prototype 
of a situation where a phase is uniquely defined in an 
obvious way. 

We now ask ourselves whether all three approaches give 
in general the same phase. If this is not the case, then 
each approach might be useful by itself, but it cannot 
claim to be really fundamental and universal. 

In Fig. 2 we show part of a signal. A phase portrait 
obtained by plotting x{t) against x{t + t) is shown in 
Fig. 3, a similar phase portrait using the Hilbert trans- 
form in Fig. 4. In neither case one sees a point around 
which the orbit circles, thus neither allows a clear and 
robust definition of phase. The spectrum, obtained with 
a Welch window, is shown in Fig. 5. A prominent peak is 
seen, but this peak is not sharp and thus a unique angu- 
lar frequency seems not obtainable. One can try several 
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FIG. 3: Phase portrait of the signal shown in Fig. 2, obtained 
by plotting x{t) against x{t + r) with t = 0.015). 



other methods popular in signal analysis, but we argue 
that none of them will lead to a robust determination of 
a phase. 

And yet - there is a simple and clear-cut phase that 
enters in this example. The signal shown in Fig. 2 is gen- 
erated by a random process defined as 

x{t) = cos((/)(t)) (2) 

with the phase performing a biased random walk (cf. 

d(t){t)/dt ^ LU + r](t) (3) 
where r](t) is (5-correlated white noise, 

(77(t)) = 0, mv{t')) = D5{t - <')• (4) 

The parameter values used in Figs. 2 to 5 are w = 1 Hz 
and D = 5, and the integration was made with step 6t — 
0.0001. The delay used in Fig.3 was r = 0.015. 

If the noise variance D were much smaller, we would 
not have much problems. The problems arise since we 
chose a rather large D such that the phase is not mono- 
tonically increasing. Instead there are large intervals dur- 
ing which the phase decreases, leading to "fake" loops in 
Figs. 3 and 4. Our point is not that presently popu- 
lar methods for extracting phases cannot distinguish be- 
tween such phase reversals (or even just sudden slow- 
downs of the instantaneous phase velocity) and "true" 
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FIG. 4: Phase portrait of the signal shown in Fig. 2, obtained 
by plotting x{t) against its Hilbert transform. 
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FIG. 5: Spectrum of the signal shown in Fig. 2. 



amplitude variations |[9|. Rather, we want to stress that 
there is no way in principle to distinguish between them. 
Thus attempts to improve on phase extraction methods 
in similarly ambiguous situations ITU are likely to 
lead to ambiguous results, even if this ambiguity might 
be hidden. 

Ways out of avoiding these ambiguities can be found 
only by restricting what we accept as a sensible phase 
definition. One could argue e.g. that a basic intuitive 
feature of a phase is its continuous temporal progression, 
i.e. positivity of the instantaneous phase velocity. De- 
manding this would mean that there is no possibility at 
all to define a phase for the above model, and the same 
would be true for a large class of signals. 

Does this mean that such a requirement is too restric- 
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tive to be useful? We believe not. One traditional way 
out of the dilemma when phases should be defined for ar- 
bitrary signals is Fourier analysis. One decomposes the 
signal into harmonic components, and can then define 
phases for each component (or, when the signal is de- 
composed into frequency bands, for each band). What 
we propose is to decompose signals more generally into 
components with positive but not necessarily constant 
(as in a Fourier decomposition) phase velocities. This 
added freedom might allow much more physically rele- 
vant decompositions. Indeed, we do not have to invent 



any new example for this, since the best example demon- 
strating the power of such an approach is known since 
nearly four hundred years: progress in understanding 
planetary motion was only possible when Kepler replaced 
the decomposition into the harmonic epicycles of Ptole- 
maeus and Copernicus by a decomposition into elliptic 
motions, which are just of the type advocated by us 
Details of such a decomposition will of course depend on 
the problem at hand, and we cannot give any general al- 
gorithm. But the possibility and the eventual usefulness 
of such an approach should be kept in mind. 
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